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Abstract 



X 

Pk 

» I . Strong astrophysical shocks, diffusively accelerating cosmic rays (CR) ought to develop CR 

"£i ' precursors. The length of such precursor L p is believed to be set by the ratio of the CR mean free 

C^J , path X to the shock speed, i.e., L p ~ cX/V s i, ~ cr g /V s i,, which is formally independent of the 

CR pressure P c . However, the X-ray observations of supernova remnant shocks suggest that the 
precursor scale may be significantly shorter than L p which would question the above estimate 
unless the magnetic field is strongly amplified and the gyroradius r g is strongly reduced over 
a short (unresolved) spatial scale. We argue that while the CR pressure builds up ahead of the 
shock, the acceleration enters into a strongly nonlinear phase in which an acoustic instability, 
f*^) , driven by the CR pressure gradient, dominates other instabilities (at least in the case of low /3 

ff*j ' plasma). In this regime the precursor steepens into a strongly nonlinear front whose size scales 

r — , with the CR pressure as Lf ~ L p ■ \L S /L p ) (P c /Pg) , where L s is the scale of the developed 

f ^ ■ acoustic turbulence, and P c /P g is the ratio of CR to gas pressure. Since L s <g L p , the precursor 

f"*) | scale reduction may be strong in the case of even a moderate gas heating by the CRs through the 

acoustic and (possibly also) the other instabilities driven by the CRs. 
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1 Introduction 

Significant progress in observations of galactic supernova remnant shocks (SNR) over 
the passed decade have furthered our understanding of the particle acceleration mech- 
anism which is deemed responsible for both the individual SNR emission and for the 
galactic cosmic ray (CR) production as a whole. The mechanism is known as the first 
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order Fermi or diffusive shock acceleration (DSA, see [12] for a comprehensive and 
12611 for short review). The shocks are typically observed as surprisingly thin filaments, 
particularly well resolved in the X-ray band. The working hypothesis is that this emis- 
sion is due to the super- TeV shock accelerated, synchrotron radiating electrons, occa- 
sionally also seen in gamma-rays (through the inverse Compton (IC) up-scattering of 
the background photons). The gamma emission may in some cases be contaminated or 
even dominated by the accelerated protons via n° decay QjJ, |2|] . A convincing demon- 
stration of namely the latter scenario would, of course, be a prima facie evidence for 
the acceleration of also the main, i.e. hadronic component of galactic CRs in SNRs. 

The key element of the DSA is multiple crossing of the shock front with a ~ V,/,/c 
energy gain after each crossing. In doing so particles diffusively escape from the shock, 
on average to a distance L p ~ K (p) /V s t,. Here K is the momentum dependent diffusion 
coefficient and V s k is the shock velocity. An obvious morphological consequence of this 
process should be and extended ~ L p shock precursor filled with synchrotron radiating 
electrons. A disturbing fact of the high-resolution X-ray observations is that there is 
no supporting evidence for such extended precursors [28, atZLlal- Note that L p grows 
with particle momentum as almost certainly does K(p), so particles of higher energy 
should make thicker emission filaments than do the lower energy particles. This trend 
does not seem to be supported by the observations either, see IJ7L I14I1 . 

According to another widely accepted view, the particle diffusion coefficient K 
should be close to the Bohm value, K ~ cr g (p) /3, which requires strong magnetic 
fluctuations 8B^ ~ Bo at the resonant scale k ~ \/r g (p). The high level of fluctua- 
tions is achieved through one of the instabilities driven by accelerated particles. The 
following three CR driven instabilities have been suggested to generate magnetic field 
fluctuations. The first one is the well known [36] ion cyclotron resonant instability of a 
slightly anisotropic (in pitch angle) CR distribution ahead of the shock |9]. The free en- 
ergy source of this instability has the potential to generate very strong field fluctuations 

(8B/B ) 2 ~M A P c /pVl (1) 

where Ma 3> 1 is the Alfvenic Mach number, P c is the CR pressure, p is the gas density 
and u s is the shock velocity. However, the actual turbulence level was thought to remain 
moderate, 8B ~ B ( e.g., MH). 

The second instability, is a nonresonant instability driven by the CR current. The 
advantage of this instability seems to be twofold. First, it cannot be stabilized by the 
quasilinear deformation of the CR distribution function since in the upstream plasma 
frame the driving CR current persists, once the CR cloud is at rest in the shock frame. 
Second, it generates a broad spectrum of waves, and the longest ones were claimed 
to be stabilized only at the level 8B ^> Bq, due to the lack of an efficient stabilization 
mechanism at such scales. Within the context of the CR acceleration, this instability 
was first studied by Achterberg [4], but the potential of this instability to generate 
strong magnetic fields was first emphasized in refs. Bell & Lucek 111 ill . Bell II 1 Ofl , so 
the instability is often referred to as the Bell's instability. The stabilization mechanism 
was assumed to be the magnetic tension [10]. Bell also pointed out that in the most 
interesting regime, the instability is driven by a fixed CR return current through the 
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Ampere force J c x B. It should be noted, however, that the dissipation of the return 
current due to the anomalous resistivity still needs to be addressed. 

The third instability is an acoustic instability (also known as Drury's) driven by the 
pressure gradient of accelerated CRs upstream [20]. This instability was also studied 
numerically by Dorfi 111 811 . The pressure gradient is clearly a viable source of free 
energy for the instability. So, among quantities varying across a shock, the pressure 
jump is the most pronounced one in that sense that it does not saturate with the Mach 
number, unlike the density or velocity jumps. 

Curiously enough, the acoustic instability has received much less attention than the 
first two. Moreover, in many numerical studies of the CR shock acceleration, special 
care is taken to suppress it. The suppression is achieved by using the fact that a change 
of stability occurs at that point in the flow where <9lnK"/<9lnp ~ — 1 (for both stable 
and unstable wave propagation directions). Here K is the CR diffusion coefficient, and 
p is the gas density. Namely, one requires this condition to hold identically all across 
the shock precursor, i.e., where the CR pressure gradient VP C j£ 0. Not only is this re- 
quirement difficult to justify physically, but, more importantly, an artificial suppression 
of the instability eliminates its genuine macroscopic and microscopic consequences, as 
briefly discussed below. 

Among the macroscopic consequences an important one is the vorticity genera- 
tion through the baroclinic effect (missalignment of the density and pressure gradients 



Vp x VP ^ 0, e.g. (270). Here VP may be associated with a quasi-constant macro- 
scopic CR-gas pressure gradient VP C , generally directed along the shock normal. Vari- 
ations of Vp are locally decoupled from P c , unlike in the situation in a gas with a 
conventional equation of state where P — P(p) and where the baroclinic term van- 
ishes. The vorticity generation obviously results (just through the frozen in condition) 
in magnetic field generation, so that the field can be amplified by the CR pressure gra- 
dient. More importantly, this process amplifies the large scale field, required for accel- 
eration of high energy particles. Furthermore, the amplification takes place well ahead 
of the gaseous subshock. The both requirements are crucial for improving high energy 
particle confinement and making the shock precursor shorter in agreement with the ob- 
servations. Large scales should be present in the ambient plasma as a seed for their 
amplification by the acoustic instability and could be driven/seeded by wave packet 
modulations. Apart from that, they result from the coalescence of shocks formed in the 
instability, and from the scattering of Alfven waves in fc-space by these shocks to larger 
scales [17]. Note that the Bell instability is essentially a short scale instability (the 
maximum growth rate is at scales much smaller than the gyro-radii of accelerated par- 
ticles). At larger scales the growth rate decreases and the instability transitions into the 



conventional cyclotron instability. It should be noted that, as well known 1251127112111 . 
vorticity (and thus magnetic field) can be efficiently generated also at the subshock. 
This would, of course, be too late for improving particle confinement and reducing the 
scale of the shock precursor. 

Now the question is which instability dominates the CR dynamics? Given the finite 
precursor crossing time, it is reasonable to choose the fastest growing mode and con- 
sider the development of a slower one on the background created by the fast mode after 
its saturation. The resonant cyclotron instability is likely to dominate at the outskirt of 
the shock precursor where both the CR current and pressure gradient (driving the other 
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two instabilities) are weak, whereas the pitch angle anisotropy is strong enough to drive 
the resonant instability. Recall that the anisotropy is typically inversely proportional to 
the local turbulence level which must decrease with the distance from the shock. In 
this paper, however, we focus on the main part of the shock precursor where both 
the CR-pressure gradient and CR current are strong, and consider the Bell and Drury 
instabilities as the strongest candidates to govern the shock structure. In fact, these 
instabilities are coupled, not only by the common energy source but also dynamically. 
The coupling of the magnetic and the density perturbations, akin to the modulational 
instability will be the subject of a separate publication. This paper will be limited to 
a simpler situation in which one of the instabilities dominates. In particular we will 
identify conditions under which the acoustic instability grows more rapidly. Then, we 
determine the shock structure resulting from the nonlinear development of the acoustic 
instability. 

2 Comparison of the growth rates 

The derivations of the linear dispersion relations of the Bell and Drury modes can be 
found in refs.[20, lOfl. For the magnetic field and density perturbations of the form, B <* 
exp(—ikx ± icot) and p °= exp (ikx ± icot), in an approximate symmetric representation, 
in which the diffusive CR damping of acoustic mode and the resonant CR contribution 
to the nonresonant magnetic instability are neglected, these relations can be written as 
follows 



co 2 = k 2 C 2 -2y B kC A (2) 

co 2 =k 2 C 2 -2iy D kC s (3) 

Here Ca and C$ are the Alfven and the sound speeds, respectively while 

P* , 1 dP c 

V Po 2p C 5 dx 

The last expressions approximately represent the maximum growth rates of the non- 
resonant B-field (Bell) perturbations and acoustic density (Drury) perturbations, driven 
by the CR return current and by the CR pressure gradient, respectively. Note that the 
instabilities are different in that the Bell mode is unstable in the limited A:-band where 
< k < 2/b/Ca, and Yb is only the maximum growth rate that is obviously achieved at 
k = Yb/Ca, whereas the acoustic mode saturates with k at k > Yd/Cs at the level Yd- 

Next, we compare the growth rates of the two instabilities. The comparison should 
be done for the same equilibrium distribution of CRs upstream of the subshock, i.e., in 
the CR precursor. Apparently, there is a difficulty here. The acoustic instability analy- 
sis presumes a certain level of magnetic fluctuations (e.g. 8B ~ B) to pondermotively 
[3] couple CRs to the gas flow, and thus ensure the equilibrium. In the Bell's stabil- 
ity analysis, it is assumed that the undisturbed B-field is along the shock normal and 
Jcr points into the same direction. Then, the CR pressure gradient cannot be statically 
compensated since Jcr xB = 0. Implicitly though, one may assume that the resonant 
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instability provides Alfven waves which balance the CR pressure gradient along the 
main field through the wave pondermotive pressure at larger scales. Then, the Bell in- 
stability develops at the scales much shorter than the gyro-radii of the current-carrying 
particles. Note that such treatment cannot be extended to the larger scales without fur- 
ther considerations B17lll3ll . In a strongly modified shock, the equilibrium distribution 
upstream in a steady state is given by [30]: 



/ = /o(p)exp 



q{p) 



0< 
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(5) 
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where is the flow potential, u = d(j)/dx < 0, /o (p) is the CR distribution and q (p) = 
—dlnfo/dlnp the spectral index at the subshock. On the other hand, if the shock 
modification is negligible, the equilibrium is simply 



/ = /o 0) exp 



1 



K(p) 



<j>(x) 



x>0 



(6) 



One sees that the only difference between the last two representations of the parti- 
cle distribution upstream is the q/3 factor in the exponent of eq.(f5J. This quantity is 
well constrained in the nonlinear solution given by eq.Q: 3.5 < q(p) < q. m b, where 
Isub = 3r«/i/ { r sub ~ 1) w i m r sub being the subshock compression ratio. For the ratio 
of acoustic to magnetic growth rates, eq.Q, we obtain 



Yd _ Ca c 2 I p 1 q 
Yb ~ C s 3co ci \ y/l+ p 2 3k- 



(7) 



where p is in units of mc. We have introduced the spectrum averaged quantity above as 
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/ ( • ) /o (/?) P 2 exp (# /3 K-) c//? 



Jfo{p)p 2 sxp(q^)/3K)dp 

and where for the test particle solution, given by eq.© one should replace q/3 
For the Bohm diffusion coefficient K = r g (p)c/3, we obtain 



(8) 
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YD = /q\CA 

Yb \3/Q 

From the last formula we may conclude that the acoustic instability dominates the 
magnetic one in the case of low j3 = C^/C\ <C 1 upstream, regardless of the degree of 
nonlinearity of acceleration. The reason for such a counter-intuitive relation between 
the growth rates of these two instabilities is that a stronger magnetic field (low-beta 
plasma) supports a stronger CR pressure gradient which drives the acoustic instabil- 
ity. Note, that plasma heating upstream would increase the role of magnetic instability 
unless the large scale magnetic field is also strengthened, e.g. through an inverse cas- 
cade of the turbulent magnetic energy [17]. On the other hand, the development of 
the acoustic instability makes (Sec|Ui the precursor shorter. This boosts the gradient 
driven acoustic instability and leaves less room for the current driven magnetic insta- 
bility. Next, we consider the CR transport in a developed acoustic turbulence with an 
admixture of Alfven waves, presumably generated by cyclotron instability at a distant 
part of the shock precursor and convected into its core. 
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We assume that magnetic perturbations in the precursor are of the following two types. 
First, there are conventional shear Alfven waves, stemming from the CR cyclotron 
instability. Second, there are compressible magnetosonic perturbations generated by 
the Drury instability. The pitch angle scattering of CRs in the both wave fields has 
been calculated in many publications. We can use the expressions given by eqs.(6,7) 
in ref. ll 31 . After summing the Bessel function series and retaining only the magnetic 
parts of the scattering wave fields, we obtain for the pitch-angle diffusion coefficients 
the following expressions 



K = -(l-M 2 )E^// A (^x,r) 

x cos {k\\c pi i ) a i ^— ^ J $ 2<;sin— - 

11 ot z V 2 
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Here I A and I s are the normalized (to Bq/%%) spectral densities of magnetic fluctuations 
of Alfven and slow magnetosonic wave components of the upstream turbulence, jJ, is 
the cosine of particle pitch angle, § = k±cy/l — ji 2 /Q., Jq is the Bessel function, and 
Q. is the relativistic gyrofrequency. The spatial diffusion coefficient can be evaluated 
as follows (e.g. [24]) 



K= — 



1-H< 



(12) 



where ^ = (D^ +DVj / (l — IJ- 2 ) and (•) denotes here the pitch angle averagin g. F or 
the turbulence spectra extended in k±^ (such as the Goldreich-Shridhar cascade 12211 ). 
the scattering frequency @n peaks at |ju| = 0, 1 11511 . Indeed, the scattering is known 
to be strongly suppressed for 4 3> 1 (high frequency perturbation of particle orbits), so 
that particles with |ju | ~ 1 are subjected to a more coherent (not so rapidly oscillating) 
wave field. Furthermore, particles with fi^O are effectively mirrored by the compress- 
ible component of magnetic turbulence. Once there are peaks at |ju| =0,1, ^ should 
also have minima in between, i.e. at some |/x| = jUq, where < jUn < 1. Upon writing 
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and evaluating the integral in pi by using the steepest descent method, we obtain 

__ 2%c 2 {\-^) 

^(Mo)^(Mo) 

where the double prime denotes the second derivative. Assuming that the compressible 
part of the turbulence, which originates from the acoustic instability, dominates (D 4 <C 
D s ), we introduce the the following momentum averaged diffusion coefficient, that will 
be used in the next section 

_ i= fK(p)f(p)p 3 d P = Kg 
If(p)p 3 d P F s + a 

Here Kb is the momentum averaged Bohm diffusion coefficient Kb = cr g (p) /3, a is 
the normalized level of Alfvenic turbulence a ~ {8Ba) jB\, (originating from Efi), 
and Fs is the level of compressible turbulence (originating from D s ). The latter can be 
calculated in the simplest ID model in which the acoustic waves unstably grow at a 
rate Jd and then steepen into an ensemble of shocks (shocktrain) H3lL 13311 



dPcV . ( U x -' 



Fs = <*{-£) with o=[j^ s ) &i»n B ) (14) 

Here L s is the average distance between the shocks in the ensemble, and & ~ 1 depends 
predominantly on the shock obliquity and on the specific model for the shocktrain 
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4 Turbulent CR Front Solution 

In this section we derive the equations for the shock transition from the kinetic convection- 
diffusion equation that can be written in a steady state as follows 

d_f dg\ _]_du dg_ 

dx\ dx J 3 dx dp 

where g = fp 3 and u is the (positive) magnitude of the flow speed. We assume that 
the CR pressure integral is dominated by particles with p ^> mc, but there is a spectral 
break at some p = p* which makes the integral finite without introducing a cut-off 
momentum 13 ill . This makes the well known two-fluid model [19] suitable for our 
purposes. More problematic is the introduction of Kb, eq.([T3"l>. which may require 
the upper cutoff. Integrating the appropriately weighted convection-diffusion equation 
( TT3T > in momentum, and supplementing the result with the conservation of the mass and 
momentum fluxes across the shock transition, we arrive at the following closed system 
of equations that governs the shock transition 

d ( k b dP c \ 1 du 

^-\UP C + =-5 ^ = -- -^-P c 

dx\ a {dP c /dx) z + a ox J 3dx 
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pu = p\U\ = const 



pu + P C = pi«! = const (16) 

Here P c is the CR pressure, p\ and u\ — V S h are the far upstream density and flow speed. 
From these equation we obtain the following equation for the shock front structure 

dP c \ 2 , L B dP c , a 



dx ) oP c (l-P c /P c2 ) dx a 

where P& = {6/l)p\u\, Lg = Kb/u\ and we assumed that thermal gas is not heated 
appreciably, so its pressure can be neglected in eq.(fT6t [30]. The shock solution can 
be obtained from the last equation in a closed form as x = x (P c ). Introducing a dimen- 
sionless coordinate z = xLg/oP% an d 3* = ?-P c /P C 2 — 1 we can rewrite the last equation 
describing the front transition as 

<94>\ 2 8 .90 

4a = 



dz J 1 - 4> 2 dz 

where the transition is governed by a single parameter a = a<3 (Pc-i/Lb) ■ Obviously, 
a smooth transition exists only for 0<a<4. Ifa^Cl<I> = — tanh (az/2), so that the 
scale of the front is the familiar Lg — Kfc/»i. If, however, a ~ 1, the scale of the shock 
(front) transition reduces to 

Lf -4L B ^ 2 -4L B {pC^) ■ 

5 Conclusions 

We have obtained a one parameter family of smooth, strongly nonlinear shock front 
transitions which are typically significantly shorter than the conventional (~ k/Vj/,) 
shock precursors. The phenomenon in general is somewhat similar to the transport 
bifurcation phenomenon, which is a subject of active research in magnetic fusion (L-H 
bifurcations, lE3l.fl6l.l32h. 



The critical parameter that governs such transitions (parameter a, SecJUi depends 
(through the growth rate of the front supporting acoustic turbulence, Yd) on the thermal 
gas pressure inside the front, de facto on the turbulent heating efficiency. The smooth 
transition exists for a < 4, so that in the case a > 4, a gaseous subshock must form. 
This however, should not result in a longer precursor. The determination of the heating 
efficiency, and thus a more complete study of the shock fronts obtained in this paper 
will be the subject of a separate publication. 
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by the Department of Energy, Grant No. DE-FG02-04ER54738 is gratefully acknowl- 
edged. 



5 Conclusions 



References 



I 

[20 
[21 
[22 
[23 
[24 
[25 



Abdo, A. A., & et al. 2010, Science, 327, 1 103 

Abdo, A. A., & et al. 2010, Astrophys. J., 712, 459 

Achterberg, A. 1981, Astronomy and Astrophys., 98, 195 

— . 1983, Astronomy and Astrophys., 1 19, 274 

Achterberg, A., & Blandford, R. D. 1986, Mon. Not. R. Astron. Soc. , 218, 551 

Ballet, J. 2006, Advances in Space Research, 37, 1902 

Bamba, A., Yamazaki, R., Ueno, M., & Koyama, K. 2003, Astrophys. J., 589, 827 

Bamba, A., Yamazaki, R., Yoshida, T., Terasawa, T., & Koyama, K. 2005, Astro- 
phys. J., 621,193 

Bell, A. R. 1978, Mon. Not. R. Astron. Soc. , 182, 147 

— . 2004, Mon. Not. R. Astron. Soc. , 353, 550 

Bell, A. R., & Lucek, S. G. 2001, Mon. Not. R. Astron. Soc. , 321, 433 

Blandford, R., & Eichler, D. 1987, Phys. Rep. , 154, 1 

Bykov, A. M., Osipov, S. M., & Toptygin, I. N. 2009, Astronomy Letters, 35, 555 

Cassam-Chenai', G., Hughes, J. R, Ballet, J., & Decourchelle, A. 2007, Astrophys. 
J., 665, 315 

Chandran, B. D. G. 2000, Physical Review Letters, 85, 4656 

Diamond, P. H., Liang, Y.-M., Carreras, B. A., & Terry, P. W. 1994, Physical 
Review Letters, 72, 2565 

Diamond, P. H., & Malkov, M. A. 2007, Astrophys. J., 654, 252 

Dorfi, E. 1984, Advances in Space Research, 4, 205 

Drury, L. O., & Voelk, J. H. 1981, Astrophys. J., 248, 344 

Drury, L. O. C, & Falle, S. A. E. G 1986, Mon. Not. R. Astron. Soc. , 223, 353 

Giacalone, J., & Jokipii, J. R. 2007, Astrophys. J. Lett. , 663, L41 

Goldreich, P., & Sridhar, S. 1995, Astrophys. J., 438, 763 

Hinton, F. L. 1991, Physics of Fluids B, 3, 696 

Jokipii, J. R. 1966, Astrophys. J., 146, 480 

Kevlahan, N. 1997, Journal of Fluid Mechanics, 341, 371 



5 Conclusions 10 

[26] Kirk, J. G., & Dendy, R. O. 2001, JOURNAL OF PHYSICS G-NUCLEAR AND 
PARTICLE, 27, 1589 

[27] Kulsrud, R. M., Cen, R., Ostriker, J. P., & Ryu, D. 1997, Astrophys. J., 480, 481 

[28] Long, K. S., Reynolds, S. P., Raymond, J. C, Winkler, P. E, Dyer, K. K., & Petre, 
R. 2003, Astrophys. J., 586, 1 162 

[29] MacKenzie, J. E, & Voelk, H. J. 1982, Astronomy and Astrophys., 1 16, 191 

[30] Malkov, M. A. 1997, Astrophys. J., 485, 638 

[31] Malkov, M. A., & Diamond, P. H. 2006, Astrophys. J., 642, 244 

[32] — . 2008, Physics of Plasmas, 15, 122301 

[33] — . 2009, Astrophys. J., 692, 1571 

[34] Malkov, M. A., Kennel, C. E, Wu, C. C, Pellat, R., & Shapiro, V. D. 1991, 
Physics of Fluids B, 3, 1407 

[35] Medvedev, M. V., & Diamond, P. H. 1996, Physics of Plasmas, 3, 863 

[36] Sagdeev, R. Z., & Shafranov, V. D. 1961, Soviet Phys. JETP, 12, 130 



